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Abstract 



We show that the hard-square lattice gas with activity z = —1 has a number of remarkable 
properties. We conjecture that all the eigenvalues of the transfer matrix are roots of unity. 
They fall into groups ("strings") evenly spaced around the unit circle, which have interesting 
number-theoretic properties. For example, the partition function on an M x TV lattice with 
periodic boundary condition is identically 1 when M and N arc coprime. We provide evidence 
for these conjectures from analytical and numerical arguments. 

1 Introduction 

The hard-square model is a weh-known model of two-dimensional statistical mechanics 012]. It 
describes a classical gas of particles on the square lattice, with the restriction that particles may 
not be on adjacent sites. The activity z is the Boltzmann weight per particle. One can think of 
the particles as hard squares with area twice that of a lattice plaquette. The squares are placed 
with their centers on a lattice site and their corners at the four adjacent sites. The restriction 
amounts to not allowing the squares to overlap, although they can touch. The particles/squares 
do not interact, except via this hard core. The partition function Z for the hard square model is 
then simply 



where A{n) is the number of allowed configurations with n particles. Combinatorialists would 
describe A{n) for an M x A*" lattice as the number of M x matrices with n entries 1 and the 
remaining zero, such that no row or column contains two consecutive nonzero entries. 

Despite its simple definition, the hard-square model has a number of interesting properties 
PI. In this paper, we discuss the fascinating structure present for the special case z = —1. We 
will present substantial analytic and numerical evidence that when z = —1, the eigenvalues of 
the transfer matrix with periodic boundary conditions are all roots of unity. Moreover, for an 
M X N lattice with periodic boundary conditions in both directions, we find Z = 1 when M and 

are coprime. 




(1) 
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There are a number of reasons why it is interesting to study the hard square model at negative 
activity, even though the model has negative Boltzmann weights. 

First, gases with hard cores are generically expected to exhibit a phase transition in the 
Yang-Lee universality class at some negative value of z j^lS]. For the hard-square model, this 
transition takes place at z = Zc ~ — .1. Thus our z = —1 results are describing the regime "past" 
this transition. This regime is very poorly understood, even though it should be described by an 
integrable field theory |6j (the Yang-Lee conformal field theory with a perturbation of opposite 
sign than usual). In this regime there are level crossings as z is decreased past Zc [3, making a 
field-theory analysis difficult. 

Second, the hard-square model is not known to be integrable for any values of z except for the 
trivial cases z = and z = oo. We thus do not know the origin of the behavior discovered here, 
although the degeneracies of the levels and other results we describe below do hint that there are 
symmetries yet to be uncovered. Hopefully such symmetries will be useful in understanding the 
hard square model for values of z other than —1. 

Third, gases with negative activity have been shown rigorously to be equivalent to branched 
polymers in two dimensions higher [H], and lattice animals in one dimension higher [51. More 
precisely, the partition function of the lattice gas at negative z is the generating function for 
branched-polymer configurations. 

Fourth is that the partition function of z = — 1 lattice gases arises very naturally in the study 
of an interesting class of quantum models with supersymmetry [2]. The simplest such model 
consists of interacting fermions hopping on a lattice, subject to the constraint that they cannot 
occupy adjacent sites. The partition function of the classical hard-core model on the same lattice 
at z = — 1 is the Witten index of this quantum theory, and as such gives a lower bound on the 
number of ground states of the theory. Thus the hard-square partition function at z = — 1 is the 
Witten index for this supersymmetric model on the square lattice. We will discuss this and other 
two-dimensional supersymmetric lattice fermions in a companion paper jlOj . 

Finally, we believe that the results described in this paper alone justify the study of this 
model: we know of no other non-trivial model of two-dimensional statistical mechanics whose 
transfer matrix obey the intriguing properties described below. For example, the hard-hexagon 
model (the analogous model on the triangular lattice) is integrable, but numerical results indicate 
that for no value of z do its transfer-matrix eigenvalues become roots of unity. 

In section 2 we define the transfer matrix used to obtain these results. Readers interested in 
the results can safely skip most of this and proceed to section 3, where our analytic and numeric 
results for the partition function of the hard-square model at z = — 1 are described. 

2 Transfer Matrix 

It is convenient to study the hard-square model by using its transfer matrix. For simplicity, we 
will consider periodic boundary conditions, although much of what we say in this paper applies 
to open boundary conditions as well. 

The dimension of the transfer matrix D]\f is the number of allowed configurations along a circle 
with N sites. We index these configurations by an integer i = 1 . . . D^; each configuration i with 
p particles is specified by the p integers (ii, Z2, . . . ip), which give the locations of the particles in 
this configuration. The hard core means that ir ^ is, is i l(modA^) for any r or s. The number 



moving from one site to the next. The contribution to Djv of configurations with p particles is 



of configurations around a circle is found by diagonalizing the transfer 
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the coefHcient of if in 



This generating function obeys the recursion relation 

DN{y) = D^-iiy) + yDN_2{y) 

The total number of configurations allowed with any p is Dn = -DAr(l)- For example, Di = 1 
(coming from the configuration with no particles) and D2 = 3 (one configuration with no particles, 
and two with one particle; because of the hard core two particles cannot be on two consecutive 
lattice sites). Thus this sequence of Dn goes as 1,3,4,7,11,18,. . . ; see table 1 below. These are 
called Lucas numbers, and obey the same recursion relation as do Fibonacci numbers. If desired 
one can enumerate the configurations for any N iteratively by this sort of approach. For example, 
the number of lattice configurations with two consecutive fixed sites empty are the Fibonacci 
numbers 1,1,2,3,5,... ; the number with one of those two sites occupied is = 

Consider two configurations i, j having pi and pj particles respectively. These two configu- 
rations (each around a circle) arc allowed to be next to each other in the full two-dimensional 
model if ir ^ js for all r and s. The transfer matrix Tjv acts on a vector space C^^; each basis 
element Vi of this vector space corresponds to a configuration i, and has 1 in the iih place and 
zeroes otherwise. The partition function with periodic boundary conditions in both directions is 
then 

Z{M,N)=tT (Tn)^ . (2) 

For the hard square model, if ir = js for some r and s, then the transfer matrix entry Tij = 0. If 
the two configurations are allowed next to each other, then 

(TN)ij = i allowed next to j (3) 

If z is negative, by convention we take the positive sign of the square root. The transfer matrix 
is not unique, but will yield the same partition function for any exponent Xpi + (1 — X)pj; the 
above definition makes the matrix symmetric. 

Since the boundary conditions are periodic, the model has a translation symmetry. The 
translation generator acts on the vector space C^^ as well, taking a configuration with particles 
at (ij, 12, ■ ■ ■ ip) to the configuration {ii + l,i2 + 1, ■ ■ - ip + 1) where all locations are interpreted 
mod A^. Thus = 1, and the eigenvalues t of T obey = 1. It is easy to verify that 
[Tn,T] = 0. The transfer matrix therefore breaks into blocks T^it) acting on the eigenstates 
of T with eigenvalue t. Working in this basis makes numerical computations more tractable, 
reducing the sizes of the matrices involved by roughly a factor of N. (Eventually this doesn't 
help much, because grows exponentially.) 

An eigenstate of T with eigenvalue t can be formed from each Vi via 

V[i]{t) = {vi + t-^Tvi + t-^r\ + ■■■ + i-(^-i)r^-ii;,) V^/AT, 

where A/^ is the smallest integer which has T^^Vi = vi. The normalization is chosen so that 
Vj*j • V[j] = 1. For example, the state with no particles has Mi = 1, while the state i with particles 

at (2, 4, . . . A^) for even has Mi = 2. The state V[j](t) is non-zero only if t^' = 1. Even though 
the transfer matrix does not conserve the number of particles, T does, so the states V[j](t) have 
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a fixed number of particles pi. Obviously any state Vj which obeys Vj = T^Vi for some integer 
r results in Vy] = f'V[i]. Thus to give a complete (but not over-complete) set of states we must 
choose just one particular i in each of these sets. 

Let us first examine the action of the transfer matrix in the sectors with t ^ 1. The simplest 
state in these sectors is the one-particle state, which we denote as [1] (the translation-invariant 
state with no particles is nonzero only in the t = 1 sector). The transfer matrix takes a onc- 
particle state with a particle on the rth site to a linear combination of all states which do not 
have a particle in the (r)th place. Using this, we see that the matrix element 

iTN)[l][l]=zit + f + ...t'') = -Z, 

where we used "^^=0^-^ ~ ^- general matrix elements, we need the function T([i];t) = 
Ylr=i'^^'^ ^'^^ ^^(^^ state [i], where as before the configuration i has particles at (ii, ^2, . . . , ^p). 
Then 

(T^.)[,][i](t^l) = ->+i)/2y^r([z];t) 

while (T/v)[i][j](t) = (T/v)[,j][]^](l/t). Each term in these matrix elements arises when configurations 
are forbidden to be next to each other. To get the general matrix elements, the idea is likewise to 
see which configurations are forbidden. The end result is related to the product T([i]; t)T([j]; l/t), 
but to not overcount forbidden configurations, each term must have coefficient 1. Precisely, by 
using = 1 rewrite the product as 

N-l 



k=0 

Then we have 

{T^)^M(t ^ 1) = -z^^^p^y^^^^ J2 (4) 

where 9{a) = 1 — 6ao, i-e. ^(0) = and is 1 otherwise. Note that if z is real and positive, Tjy is 
Hermitean. In the case of interest here z is not, but if desired one can redefine TN{t) without 
changing its eigenvalues to make it symmetric. 
Similar arguments give Tjv(t = 1): 



(T^)[,]L,](t = 1) = (n - X;^(a,([z], bD)) (5) 

\ k=0 / 



where the are defined from T([i]; i)r([j]; l/t) as above. Note that when [i] is the configuration 
with no particles, T([i]; t) = and so all the corresponding = as well. 

3 The partition function 

We write the partition function in terms of the roots of the characteristic polynomial of the 
transfer matrix defined in the previous section. The characteristic polynomial Pn{x) is 
defined as 

Pn{x) =det{x-TN). (6) 
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N 

~r 

2 
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4 

5 
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7 
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10 

11 

12 

13 

14 

15 



Dn 
1 
3 
4 
7 
11 
18 
29 
47 
76 
123 
199 
322 
521 
843 
1364 



l)2(x 



l)4(a;22 
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X - 1 
(x2 + l)(a; - 1) 
(x^ - l)(a; - 1) 
{x^ - l){x^ - l){x - 1) 
{x^ + l)2(x - 1) 

(x^ - l)2(a;3 - l)(x2 + l)(x - 1) 

{x^^ + l)\x-l) 

(xlS-l)2(x9-l)4(x3-l)(x-l) 
(xl4 - l)5(xS - l)5(x5 - l)2(x2 + l)(x - 1) 
(x55-l)2(x22 + 1)4(3; _1) 
1)6(^12 _ 1)12(2.6 _ 1)2(^,4 _ i)(^3 _ i)(^2 



I)4(x4-l)(x2-l)(x-l) 
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l)(x-l) 



■ l)7(xi6 
1)6 (x^s . 



I)4(x2^ - l)4(xl'^ + 1 



1) 



1)28(^14 
1)18/^15 _ 



- I)2(xl0 - l)7(x7 + l)4(x2 + l)(x - 1) 
I)12(x5 + l)2(x3-l)(x-l) 



Table 1: Characteristic polynomials of the transfer matrices T/v. The roots of this polynomial 
give the partition function for an M x lattice with periodic boundary conditions, as in Q. 

The partition function for an M x A'" lattice with periodic boundary conditions in both directions 
is 

ziM,N) = ^{xmr- (7) 

i=l 

where the x, are the roots of Pn{x). Because z is negative, the transfer matrix is not Hermitean; 
a resulting complication we will discuss is that not all roots Xj need be eigenvalues of T^. 

Our main conjecture is that the roots Xi{N) of the characteristic polynomial are all roots of 
unity and for a given N can be grouped into "strings" . A string is a set of Xj evenly spaced around 
the unit circle. We find just two kinds of strings, which we denote and S~ . The former are 
roots the values Xk = e*2'^^/'^ for A; = 0, 1, ... S" — 1, and the latter are the values x^ = e*'^(2fe+i)/5 
for A: = 0, 1, ... 5 — 1. The existence of a string means that the polynomial (x"^ =F 1) divides 
P]Sf{x). Moreover, we find that all the values of S for all strings for a given A'^ share a divisor 
with A^, except for a single 1"*" string for every A^. 

The strong evidence for these conjectures comes from the numerical results presented in 
Table 1. We have checked up to A^ = 15 that these conjectures hold. Another result apparent 
from this table is that det(rjv) = 1, so Pn{0) = (—1)-^^. Additional numerical results for N <9 
suggest that the roots of the characteristic polynomial for the transfer matrix with open boundary 
conditions are also roots of unity. For simplicity we will focus here on periodic. 

The reason why the roots of PNi^) are all of unit modulus is mystery to us; we do not know 
any other lattice gases even at z = —1 which share this property. Despite substantial effort, we 
have not found a formula for P/v(x) for arbitrary A". However, one can exploit these conjectures 
to better understand the string structure. We define the strings so that only one kind of string is 
present for a given S and A^. This means in P/v(x) any occurrence of (x'^ — l)(x'^ + 1) is combined 
into (x2'^ — 1). Then we let n^{N) be the number of strings for a given A^; our convention 
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means that ng{N)ng{N) = 0. Since each polynomial Pn is of order x , we have 

s 

Another useful fact is that because ^2-nikM/s _ q ^^ij^iggg M is a multiple of 5, only strings 

with M a multiple of S contribute to Z{M, N). An 5+ string contributes 5, while an S~ string 

contributes {-lY^I^S. 

More about the string structure can be learned by exploiting modular invariance, the sym- 
metry of the model Z(M, N) = Z(N, M) under interchange of the two cycles of the torus. A 
consequence is that knowing the P/v [x) for all the N less than and equal to a given S determines 
the numbers of S'-strings n|(iV) for all N. 

Let us explain how to implement this recursive procedure. Since Ti = 1, the only string for 
A'' = 1 is a single l"*" string, as apparent from Table 1. Thus Z(M, 1) = 1, and by symmetry 
Z{1,N) = 1. Because only strings where M is a multiple of S contribute to Z{M,N), only 
1-strings contribute to Z{l,N). Since Z[1,N) = 1, we have just a single l"*" string: nf{N) = 1 
and n^(N) = for all A^. This is apparent in Table 1. Moving on to A^ = 2, we have D2 = 3 
and since there is exactly one 1-string, there must be just one 2-string. By our conjecture that 
-Piv(O) = (—1)-^^, this must be a 2+-string, giving the result in the table. Thus 

Z{2j -1,2) = 1, Z{Aj - 2, 2) = -1, Z{4j, 2) = 3, j an integer > 1. 

By symmetry, the same results hold for Z(2,A^). Only 2-strings and 1-strings contribute to 
Z(2,N), and we already know there is exactly one l"*" string for all N. The explicit expression 
for Z(2, N) then tells us that 

{4:j) = (4j ~ 2) = 1, j an integer > 1, 

with all other n2 = 0. In a similar fashion, one finds n^{3j) = 1 and n'^{4j) = 1. These are all 
apparent in table 1. 

For A^ > 5, this procedure yields only part of P/v uniquely. For A^ = 5, we know there is 
one 1^ string, and no 2,3 or 4 strings. Since = 11 and ^5(0) = —1, this means that there 
are either two 5+ strings or two 5~ strings. We do not have a general conjecture which will 
distinguish between the two possibilities, so we must compute P5{x) explicitly. As seen in table 
1, there end up being two 5^ strings. Now we can use the symmetry under interchange of M and 
N to see that n^(lOj) = 2 and 715(10^ -|- 5) = 2. Likewise, to find that there are two 6"'" strings 
for N = 6, one needs to find Pe{x) explicitly. This then results in ng"(6j) = 2. The general result 
is found by noting that for all strings with S > N, Sns{N) is always a multiple of A^. This 
means that if we were using a transfer matrix in the M-dircction instead, this contribution to the 
partition function would arise from nAr(M) = Sns{N)/N strings of length A^. Using the results 
from table 1, this procedure can be applied to find 7T,y"(28j) = (14(2j -|- 1)) = 4, ng{Wj) = 5, 
n+(18i) = 6, n+{9{2j + 1)) = 4, <o(14i) = 7, and so on. 

Unfortunately, since D]\f increases exponentially with A^, the strings with S < N make only 
a relatively small contribution to the partition function. As is apparent from the table, as A^ 
increases, the number of different types of strings with S > N increases. We have not yet seen 
any pattern to these numbers, but we arc hopeful that one may exist. The increasing degeneracies 
(multiple roots at the same x) apparent in the table as N increases are an obvious hint that there 
is some yet-undiscovered symmetry structure. 



6 





1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


14 


15 


16 


17 


18 


19 


20 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


2 


1 


-1 


1 


3 


1 


-1 


1 


3 


1 


-1 


1 


3 


1 


-1 


1 


3 


1 


-1 


1 


3 


3 


1 


1 


4 


1 


1 


4 


1 


1 


4 


1 


1 


4 


1 


1 


4 


1 


1 


4 


1 


1 


4 


1 


3 


1 


7 


1 


3 


1 


7 


1 


3 


1 


7 


1 


3 


1 


7 


1 


3 


1 


7 


5 


1 


1 


1 


1 


-9 


1 


1 


1 


1 


11 


1 


1 


1 


1 


-9 


1 


1 


1 


1 


11 


6 


1 


-1 


4 


3 


1 


14 


1 


3 


4 


-1 


1 


18 


1 


-1 


4 


3 


1 


14 


1 


3 


7 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


-27 


1 


1 


1 


1 


1 


1 


8 


1 


3 


1 


7 


1 


3 


1 


7 


1 


43 


1 


7 


1 


3 


1 


7 


1 


3 


1 


47 


9 


1 


1 


4 


1 


1 


4 


1 


1 


40 


1 


1 


4 


1 


1 


4 


1 


1 


76 


1 


1 


10 


1 


-1 


1 


3 


11 


-1 


1 


43 


1 


9 


1 


3 


1 


69 


11 


A O 

43 


1 


-1 


1 


13 


11 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


12 


1 


3 


4 


7 


1 


18 


1 


7 


4 


3 


1 


166 


1 


3 


4 


7 


1 


126 


1 


7 


13 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


-51 


1 


1 


1 


1 


1 


1 


1 


14 


1 


-1 


1 


3 


1 


-1 


-27 


3 


1 


69 


1 


3 


1 


55 


1 


451 


1 


-1 


1 


73 


15 


1 


1 


4 


1 


-9 


4 


1 


1 


4 


11 


1 


4 


1 


1 


174 


1 


1 


4 


1 


11 



Table 2: The partition function Z{M,N) = Z{N,M) for M <15,N < 20. 

The conjecture that all values of S for a given N (except for a single 1^ string) share a divisor 
with S means that 

Z{M, N) = 1 when M and N are coprime. 

This fact is apparent in Table 2. For values of M and N which are not coprime, the partition 
function grows with increasing M and at a much smaller rate as generic statistical mechanical 
systems. Since all the Xi{N) have magnitude 1, the maximal value of Z(M, N) is the smaller of Z?Ar 
or Dm- D]\j grows exponentially in A^, while Z{M, N) for generic values of z grows exponentially 
in NM. Exponential growth in NM is of course the standard behavior for statistical-mechanical 
systems: it is the statement that the free energy is extensive (a notable exception is systems with 
supersymmetry) . In fact even for the analogous systems at z = — 1, the partition function grows 
exponentially with the volume. For example, for the same model on the triangular lattice (the 
hard-hexagon model), our numerics indicate that the partition function grows as ~ (1.14)^^^. 

In this section we have been careful to refer to the Xi{N) as the roots of the characteristic 
polynomial Pi\f{x), not as the eigenvalues of Tjy. The reason is that for > 4, some of these 
roots do not correspond to eigenvalues. The eigenvectors of a non-Hermitean transfer matrix 
(even one of determinant 1 like T/v) need not span the space of states if two of the roots coincide. 

'1 



A simple example is the matrix A{b) ~ yQ -^j ■ Both roots are 1, but for any 6 7^ it has 

only one eigenvector ei = . The vector cq = is linearly independent of vq, but obeys 

A{b)eo = eo + (1 + b)ei. Acting repeatedly with A does not change the coefficient of cq, but just 
continues to change the coefficient of the eigenvector ei. This is because matrices A{b) and A{b') 
commute; the associated conserved quantity is the coefficient of eo- 

Despite this (interesting) complication, the partition function for periodic boundary condi- 
tions in both directions can still be written in the familiar form ((71); the sum is over all roots, not 
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just the eigenvalues. For example, for N = 5 and t = 1, the transfer matrix is 

/ 1 i\/5 V5\ 
T^it = 1)= -4 3i 

\V5 3i 2 / 

where the states are in order of increasing particle number pi = 0,1,2. The characteristic polyno- 
mial of this is (x— l)(x-|-l)^. There is only one eigenstate of eigenvalue —1; the vector Vq = (0, 1, i) 
is orthogonal to the eigenvectors but is not an eigenvector itself. Because the coefficient of Vq is 
conserved as one acts with T5, the partition function with periodic boundary conditions remains 
(|7I). This means that in this sector the partition function is Z{M, 5)(t = 1) = — 1 for M odd, and 
Z{M,5){t = 1) = 3 for M even. However, the matrix elements of [T^(t = 1)]^ are not periodic 
in M like this. The partition function on a cylinder with fixed boundary conditions on the ends 
will not be periodic; its magnitude will continue to increase, roughly linearly in M. 

Moreover, one can use this to find a family of matrices with the same roots Xi as T5{t = 1). 
If all the roots —1,-1,1 were eigenvalues, the matrix [T^it = 1)]^ would equal the identity. 
Here Q = [T^{t = 1)]^ — 1 is non-zero; it annihilates the eigenvectors and obeys = 0. The 
characteristic polynomial of Tr,[t = 1) — XQ is independent of A. When A = 1/2, this matrix has 
three eigenvectors instead of two. Analogous results follow for all > 4. It is not clear, however, 
if the presence of this nilpotent symmetry operator Q will help in the analysis of the model. 

This work was supported in part by the foundations FOM and NWO of the Netherlands, by 
NSF grant DMR-0104799, and by DOE under grant DEFG02-97ER41027. 
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